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A METHOD AND DEVICE FOR ADAPTIVE BANDWIDTH 



PITCH SEARCH IN CODING WIDEBAND SIGNALS 




BACKGROUND OF THE INVENTION 



1. Field of the invention: 

The present invention relates to an efficient technique for 
digitally encoding a wideband signal, in particular but not exclusively a 
speech signal, in view of transmitting, or storing, and synthesizing this 
wideband sound signal. More specifically, this invention deals with an 
improved pitch search device and method. 

2. Brief description of the prior art: 

The demand for efficient digital wideband speech/audio 
encoding techniques with a good subjective quality/bit rate trade-off is 
increasing for numerous applications such as audio/video 
teleconferencing, multimedia, and wireless applications, as well as 
Internet and packet network applications. Until recently, telephone 
bandwidths filtered in the range 200-3400 Hz were mainly used in speech 
coding applications. However, there is an increasing demand for 
wideband speech applications in order to increase the intelligibility and 
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naturalness of the speech signals. A bandwidth in the range 50-7000 Hz 
was found sufficient for delivering a face-to-face speech quality. For 
audio signals, this range gives an acceptable audio quality, but still lower 
than the CD quality which operates on the range 20-20000 Hz. 

5 A speech encoder converts a speech signal into a digital 

bitstream which is transmitted over a communication channel (or stored 
in a storage medium). The speech signal is digitized (sampled and 
quantized with usually 16-bits per sample) and the speech encoder has 
the role of representing these digital samples with a smaller number of 

10 bits while maintaining a good subjective speech quality. The speech 
decoder or synthesizer operates on the transmitted or stored bit stream 
and converts it back to a sound signal. 

One of the best prior art techniques capable of achieving a 
15 good quality/bit rate trade-off is the so-called Code Excited Linear 
Prediction (CELP) technique. According to this technique, the sampled 
speech signal is processed in successive blocks of L samples usually 
called frames where L is some predetermined number (corresponding to 
10-30 ms of speech). In CELP, a linear prediction (LP) filter is computed 
20 and transmitted every frame. The L-sample frame is then divided into 
smaller blocks called subframes of size N samples, where L=kN and k is 
the number of subframes in a frame (N usually corresponds to 4-10 ms 
of speech). An excitation signal is determined in each subframe, which 
usually consists of two components: one from the past excitation (also 
25 called pitch contribution or adaptive codebook) and the other from an 
innovation codebook (also called fixed codebook). This excitation signal 
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is transmitted and used at the decoder as the input of the LP synthesis 
filter in order to obtain the synthesized speech. 

An innovation codebook in the CELP context, is an indexed set 
of A/-sample-long sequences which will be referred to as A/-dimensional 
5 codevectors. Each codebook sequence is indexed by an integer k 
ranging from 1 to M where M represents the size of the codebook often 
expressed as a number of bits b, where M=2 b . 

To synthesize speech according to the CELP technique, each 
1 0 block of N samples is synthesized by filtering an appropriate codevector from 
a codebook through time varying filters modeling the spectral characteristics 
of the speech signal. At the encoder end, the synthetic output is computed 
for all, or a subset, of the codevectors from the codebook (codebook search). 
The retained codevector is the one producing the synthetic output closest to 
15 the original speech signal according to a perceptually weighted distortion 
measure. This perceptual weighting is performed using a so-called 
perceptual weighting filter, which is usually derived from the LP filter. 

The CELP model has been very successful in encoding telephone 
20 band sound signals, and several CELP-based standards exist in a wide 
range of applications, especially in digital cellular applications. . In the 
telephone band, the sound signal is band-limited to 200-3400 Hz and 
sampled at 8000 samples/sec. In wideband speech/audio applications, the 
sound signal is band-limited to 50-7000 Hz and sampled at 16000 
25 samples/sec. 
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Some difficulties arise when applying the telephone-band 
optimized CELP model to wideband signals, and additional features need to 
be added to the model in order to obtain high quality wideband signals. 
Wideband signals exhibit a much wider dynamic range compared to 
telephone-band signals, which results in precision problems when a fixed- 
point implementation of the algorithm is required (which is essential in 
wireless applications). Further, the CELP model will often spend most of its 
encoding bits on the low-frequency region, which usually has higher energy 
contents, resulting in a low-pass output signal. To overcome this problem, 
the perceptual weighting filter has to be modified in order to suit wideband 
signals, and pre-emphasis techniques which boost the high frequency 
regions become important to reduce the dynamic range, yielding a simpler 
fixed-point implementation, and to ensure a better encoding of the higher 
frequency contents of the signal. Further, the pitch contents in the spectrum 
of voiced segments in wideband signals do not extend over the whole 
spectrum range, and the amount of voicing shows more variation compared 
to narrow-band signals. Therefore, in case of wideband signals, existing 
pitch search structures are not very efficient. Thus, it is important to improve 
the closed-loop pitch analysis to better accommodate the variations in the 
voicing level. 



OBJECTS OF THE INVENTION 



25 



An object of the present invention is therefore to provide a method 
and device for efficiently encoding wideband {7000 Hz) sound signals using 
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CELP-type encoding techniques, using improved pitch analysis in order to 
obtain high a quality reconstructed sound signal. 

SUMMARY OF THE INVENTION 

5 

More specifically, in accordance with the present invention, there 
is provided a method for selecting an optimal set of pitch codebook 
parameters associated to a signal path, from at least two signal paths, 
having the lowest calculated pitch prediction error. The pitch prediction error 

10 is calculated in response to a pitch codevector from a pitch codebook search 
device. In at least one of the two signal paths, the pitch prediction error is 
filtered before supplying the pitch codevector for calculation of said pitch 
prediction error of said one path. Finally, the pitch prediction errors 
calculated in said at least two signal paths are compared, the signal path 

1 5 having the lowest calculated pitch prediction error is chosen, and the set of 
pitch codebook parameters associated to the choosen signal path are 
selected. 

The pitch analysis device of the invention, for producing an 
20 optimal set of pitch codebook parameters, comprises: 

a) at least two signal paths associated to respective sets of pitch 
codebook parameters, wherein: 

i) each signal path comprises a pitch prediction error 
calculating device for calculating a pitch prediction error of a 

25 pitch codevector from a pitch codebook search device; and 

ii) at least one of the two paths comprises a filter for 
filtering the pitch codevector before supplying the pitch 
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codevector to the path's pitch prediction error calculating 
device; and 

b) a selector for comparing the pitch prediction errors calculated 
in the signal paths, for choosing the signal path having the lowest calculated 
pitch prediction error, and for selecting the set of pitch codebook parameters 
5 associated to the choosen signal path. 

The new method and device which achieve efficient modeling of 
the harmonic structure of the speech spectrum uses several forms of low 
10 pass filters applied to the past excitation and the one yielding higher 
prediction gain is selected. When subsample pitch resolution is used, the 
low pass filters can be incorporated into the interpolation filters used to obtain 
the higher pitch resolution. 

15 In a preferred embodiment of the invention, each pitch prediction 

error calculating device of the pitch analysis device described above 
comprises: 

a) a convolution unit for convolving the pitch codevector with a 
weighted synthesis filter impulse response signal and therefore calculating 

20 a convolved pitch codevector; 

b) a pitch gain calculator for calculating a pitch gain in response 
to the convolved pitch codevector and a pitch search target vector; 

c) an amplifier for multiplying the convolved pitch codevector by 
the pitch gain to thereby produce an amplified convolved pitch codevector; 

25 and 
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d) a combiner circuit for combining the amplified convolved pitch 
codevector with the pitch search target vector to thereby produce the pitch 
prediction error. 

In another preferred embodiment of the invention, the pitch gain 
5 calculator comprises a means for calculating said pitch gain b (J) using the 
relation: 

b V w /|y w f where ; = 0, 1, 2 7 K, and K corresponds to a 

number of signal paths, 



10 



and where x is said pitch search target vector, and y yj is said convolved 
pitch codevector. 



The present invnetion further relates to an encoder, having the 
1 5 pitch analysis device described above, for encoding a wideband input signal 
and comprising: 

a) a linear prediction synthesis filter calculator responsive to the 
wideband signal for producing linear prediction synthesis filter coefficients; 

b) a perceptual weighting filter, responsive to the wideband signal 
20 and the linear prediction synthesis filter coefficients, for producing a 

perceptually weighted signal; 

c) an impulse response generator responsive to the linear 
prediction synthesis filter coefficients for producing a weighted synthesis filter 
impulse response signal; 

25 d) a pitch search unit for producing pitch codebook parameters, 

comprising: 
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i) a pitch codebook search device responsive to the 
perceptually weighted signal and the linear prediction synthesis 
filter coefficients for producing the pitch codevector and an 
innovative search target vector; and 

ii) the pitch analysis device responsive to the pitch codevector 
5 for selecting, from the sets of pitch codebook parameters, the set 

of pitch codebook parameters associated to the path having the 
lowest calculated pitch prediction error; 

d) an innovative codebook search device, responsive to the 
weighted synthesis filter impulse response signal, and the innovative search 

1 0 target vector, for producing innovative codebook parameters; and 

e) a signal forming device for producing an encoded wideband 
signal comprising the set of pitch codebook parameters associated to the 
path having the lowest pitch prediction error, the innovative codebook 
parameters, and the linear prediction synthesis filter coefficients. 

15 

The present invention still further relates to a cellular 
communication system, a cellular mobile transmitter/receiver unit, a cellular 
network element, and a bidirectional wireless communication sub-system 
comprising the above described decoder. 

20 

The objects, advantages and other features of the present invention will 
become more apparent upon reading of the following non restrictive 
description of a preferred embodiment thereof, given by way of example only 
with reference to the accompanying drawings. 

25 
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BRIEF DESCRIPTION OF THE DRAWINGS 

In the appended drawings: 

5 Rgurel is a schematic block diagram of a preferred embodiment of 

wideband encoding device; 

Figure 2 js a schematic block diagram of a preferred embodiment of 
wideband decoding device; 

10 

JEj gure 3 is a schematic block diagram of a preferred embodiment of 
pitch analysis device; and 

^igure4Js a simplified, schematic block diagram of a cellular 
1 5 communication system in which the wideband encoding device of Figure 1 
and the wideband decoding device of Figure 2 can be used. 



20 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

As well known to those of ordinary skill in the art, a cellular 
communication system such as 401 (see Figure 4) provides a 
25 telecommunication service over a large geographic area by dividing that 
large geographic area into a number C of smaller cells. The C smaller cells 
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are serviced by respective cellular base stations 402, , 402 2 ... 402 c to 
provide each cell with radio signalling, audio and data channels. 

Radio signalling channels are used to page mobile radiotelephones 
(mobile transmitter/receiver units) such as 403 within the limits of the 
coverage area (cell) of the cellular base station 402, and to place calls to 
other radiotelephones 403 located either inside or outside the base station's 
cell or to another network such as the Public Switched Telephone Network 
(PSTN) 404. 

Once a radiotelephone 403 has successfully placed or received a 
call, an audio or data channel is established between this radiotelephone 
403 and the cellular base station 402 corresponding to the cell in which the 
radiotelephone 403 is situated, and communication between the base 
station 402 and radiotelephone 403 is conducted over that audio or data 
channel. The radiotelephone 403 may also receive control or timing 
information over a signalling channel while a call is in progress. 

if a radiotelephone 403 leaves a cell and enters another adjacent cell 
while a call is in progress, the radiotelephone 403 hands over the call to an 
available audio or data channel of the new cell base station 402. If a 
radiotelephone 403 leaves a cell and enters another adjacent cell while no 
call is in progress, the radiotelephone 403 sends a control message over the 
signalling channel to log into the base station 402 of the new cell. In this 
manner mobile communication over a wide geographical area is possible. 

The cellular communication system 401 further comprises a control 
terminal 405 to control communication between the cellular base stations 
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402 and the PSTN 404, for example during a communication between a 
radiotelephone 403 and the PSTN 404, or between a radiotelephone 403 
located in a first ceil and a radiotelephone 403 situated in a second cell. 

Of course, a bidirectional wireless radio communication subsystem 
5 is required to establish an audio or data channel between a base station 402 
of one cell and a radiotelephone 403 located in that cell. As illustrated in 
very simplified form in Figure 4, such a bidirectional wireless radio 
communication subsystem typically comprises in the radiotelephone 403: 

- a transmitter 406 including: 

10 - an encoder 407 for encoding the voice signal; and 

- a transmission circuit 408 for transmitting the encoded voice 
signal from the encoder 407 through an antenna such as 409; 
and 

- a receiver 410 including: 

15 - a receiving circuit 41 1 for receiving a transmitted encoded voice 

signal usually through the same antenna 409; and 

- a decoder 412 for decoding the received encoded voice signal 
from the receiving circuit 411. 

20 The radiotelephone further comprises other conventional 

radiotelephone circuits 41 3 to which the encoder 407 and decoder 412 are 
connected and for processing signals therefrom, which circuits 413 are well 
known to those of ordinary skill in the art and, accordingly, will not be further 
described in the present specification. 

25 

Also, such a bidirectional wireless radio communication subsystem 
typically comprises in the base station 402: 
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- a transmitter 414 including: 

- an encoder 415 for encoding the voice signal; and 

- a transmission circuit 416 for transmitting the encoded voice 
signal from the encoder 415 through an antenna such as 417; 
and 

5 - a receiver 418 including: 

- a receiving circuit 419 for receiving a transmitted encoded voice 
signal through the same antenna 41 7 or through another antenna 
(not shown); and 

- a decoder 420 for decoding the received encoded voice signal 
1 0 from the receiving circuit 419. 

The base station 402 further comprises, typically, a base station 
controller 421, along with its associated database 422, for controlling 
communication between the control terminal 405 and the transmitter 414 and 
15 receiver 418. 

As well known to those of ordinary skill in the art, voice encoding is 
required in order to reduce the bandwidth necessary to transmit sound 
signal, for example voice signal such as speech, across the bidirectional 
20 wireless radio communication subsystem, i.e., between a radiotelephone 
403 and a base station 402. 

LP voice encoders (such as 415 and 407) typically operating at 13 
kbits/second and below such as Code-Excited Linear Prediction (CELP) 
25 encoders typically use a LP synthesis filter to model the short-term spectra! 
envelope of the voice signal. The LP information is transmitted, typically, 
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every 1 0 or 20 ms to the decoder (such 420 and 412) and is extracted at the 
decoder end. 

The novel techniques disclosed in the present specification may apply 
to different LP-based coding systems. However, a CELP-type coding 
5 system is used in the preferred embodiment for the purpose of presenting a 
non-limitative illustration of these techniques. In the same manner, such 
techniques can be used with sound signals other than voice and speech as 
well with other types of wideband signals. 

10 Figure 1 shows a general block diagram of a CELP-type speech 

encoding device 100 modified to better accommodate wideband signals. 

The sampled input speech signal 1 14 is divided into successive L- 
sample blocks called "frames". In each frame, different parameters 

15 representing the speech signal in the frame are computed, encoded, and 
transmitted. LP parameters representing the LP synthesis filter are usually 
computed once every frame. The frame is further divided into smaller blocks 
of N samples (blocks of length A/), in which excitation parameters (pitch and 
innovation) are determined. In the CELP literature, these blocks of length N 

20 are called "subframes" and the A/-sample signals in the subframes are 
referred to as A/-dimensional vectors. In this preferred embodiment, the 
length N corresponds to 5 ms while the length L corresponds to 20 ms, 
which means that a frame contains four subframes (A/=80 at the sampling 
rate of 16 kHz and 64 after down-sampling to 12.8 kHz). Various N- 

25 dimensional vectors occur in the encoding procedure. A list of the vectors 
which appear in Figures 1 and 2 as well as a list of transmitted parameters 
are given herein below: 
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List of the main A/-dimensional vectors 

s Wideband signal input speech vector (after down-sampling, pre- 
processing, and preemphasfs); 
s w Weighted speech vector; 
s 0 Zero-input response of weighted synthesis filter; 
s p Down-sampled pre-processed signal; 

Oversampled synthesized speech signal; 

s T Synthesis signal before deemphasis; 

s d Deemphasized synthesis signal; 

s h Synthesis signal after deemphasis and postprocessing; 

x Target vector for pitch search; 

x' Target vector for innovation search; 

h Weighted synthesis filter impulse response; 

vy Adaptive (pitch) codebook vector at delay 7~; 

y T Filtered pitch codebook vector (v 7 convolved with h); 

c k Innovative codevector at index k (/c-th entry from the innovation 

codebook); 
c f Enhanced scaled innovation codevector; 
u Excitation signal (scaled innovation and pitch codevectors); 
u 7 Enhanced excitation; 
z Band-pass noise sequence; 
w' White noise sequence; and 
w Scaled noise sequence. 
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List of transmitted parameters 

STP Short term prediction parameters (defining A{z}); 
T Pitch lag (or pitch codebook index); 
b Pitch gain (or pitch codebook gain); 
5 j index of the low-pass filter used on the pitch codevector; 

k Codevector index (innovation codebook entry); and 
g Innovation codebook gain. 

In this preferred embodiment, the STP parameters are transmitted 
1 0 once per frame and the rest of the parameters are transmitted four times per 
frame (every subframe). 

ENCODER SIDE 



1 5 The sampled speech signal is encoded on a block by block basis by 

the encoding device 100 of Figure 1 which is broken down into eleven 
modules numbered from 101 to 111. 

The input speech is processed into the above mentioned L-sample 
20 blocks called frames. 

Referring to Figure 1, the sampled input speech signal 114 is down- 
sampled in a down-sampling module 101 . For example, the signal is down- 
sampled from 16 kHz down to 12.8 kHz, using techniques well known to 
25 those of ordinary skill in the art. Down-sampling down to another frequency 
can of course be envisaged. Down-sampling increases the coding 
efficiency, since a smaller frequency bandwidth is encoded. This also 
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reduces the algorithmic complexity since the number of samples in a frame 
is decreased. The use of down-sampling becomes significant when the bit 
rate is reduced below 16 kbit/s, although down-sampling is not essential 
above 16 kbit/s. 

After down-sampling, the 320-sample frame of 20 ms is reduced to 
256-sample frame (down-sampling ratio of 4/5). 

The input frame is then supplied to the optional pre-processing block 
102. Pre-processing block 102 may consist of a high-pass filter with a 50 Hz 
cut-off frequency. High-pass filter 102 removes the unwanted sound 
components below 50 Hz. 

The down-sampled pre-processed signal is denoted by s p (n), n=0, 1 , 
2, ...,L-1, where L is the length of the frame (256 at a sampling frequency of 
12.8 kHz). In a preferred embodiment of the preemphasis filter 103, the 
signal s p (n) is preemphasized using a filter having the following transfer 
function: 



P(z) =1 - uz~ 1 



where ^ is a preemphasis factor with a value located between 0 and 1 (a 
typical value is /x = 0.7). A higher-order filter could also be used. It should 
be pointed out that high-pass filter 102 and preemphasis filter 103 can be 
interchanged to obtain more efficient fixed-point implementations. 
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The function of the preemphasis filter 103 is to enhance the high 
frequency contents of the input signal. It also reduces the dynamic range of 
the input speech signal, which renders it more suitable for fixed-point 
implementation. Without preemphasis, LP analysis in fixed-point using 
single-precision arithmetic is difficult to implement. 

5 

Preemphasis also plays an important role in achieving a proper 
overall perceptual weighting of the quantization error, which contributes to 
improved sound quality. This will be explained in more detail herein below. 

1 0 The output of the preemphasis filter 1 03 is denoted s(n). This signal 

is used for performing LP analysis in calculator module 104. LP analysis is 
a technique well known to those of ordinary skill in the art. In this preferred 
embodiment, the autocorrelation approach is used. In the autocorrelation 
approach, the signal s(n) is first windowed using a Hamming window (having 

15 usually a length of the order of 30-40 ms). The autocorrelations are 
computed from the windowed signal, and Levinson-Durbin recursion is used 
to compute LP filter coefficients, a„ where /==1,.. M p, and where p is the LP 
order, which is typically 16 in wideband coding. The parameters a ( are the 
coefficients of the transfer function of the LP filter, which is given by the 

20 following relation: 



25 



A(z) = 1+£>,z- 1 

;=1 
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LP analysis is performed in calculator module 104, which also 
performs the quantization and interpolation of the LP filter coefficients. The 
LP filter coefficients are first transformed into another equivalent domain 
more suitable for quantization and interpolation purposes. The line spectral 
pair (LSP) and immitance spectral pair (fSP) domains are two domains in 
5 which quantization and interpolation can be efficiently performed. The 16 LP 
filter coefficients, a h can be quantized in the order of 30 to 50 bits using split 
or multi-stage quantization, or a combination thereof. The purpose of the 
interpolation is to enable updating the LP filter coefficients every subframe 
while transmitting them once every frame, which improves the encoder 
1 0 performance without increasing the bit rate. Quantization and interpolation 
of the LP fitter coefficients is believed to be otherwise well known to those of 
ordinary skill in the art and, accordingly, will not be further described in the 
present specification. 

15 The following paragraphs will describe the rest of the coding 

operations performed on a subframe basis. In the following description, the 
filter >A(z) denotes the unquantized interpolated LP filter of the subframe, and 
the filter A(z) denotes the quantized interpolated LP filter of the subframe. 

20 Perceptual Weighting: 

In analysis-by-synthesis encoders, the optimum pitch and innovation 
parameters are searched by minimizing the mean squared error between the 
input speech and synthesized speech in a perceptually weighted domain. 
25 This is equivalent to minimizing the error between the weighted input speech 
and weighted synthesis speech. 
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The weighted signal sj^ri) is computed in a perceptual weighting fitter 
105. Traditionally, the weighted signal sJin) is computed by a weighting fitter 
having a transfer function W(z) in the form: 

5 W(z)=A(z/y } ) / A(z/y 2 ) where o <y 2 <y 1 <i 

As well known to those of ordinary skill in the art, in prior art analysis-by- 
synthesis (AbS) encoders, analysis shows that the quantization error is 

1 0 weighted by a transfer function W' 1 (2) f which is the inverse of the transfer 
function of the perceptual weighting filter 105. This result is well described 
by B.S. Atal and MR. Schroeder in "Predictive coding of speech and 
subjective error criteria'*, IEEE Transaction ASSP, vol. 27, no. 3, pp. 247- 
254, June 1979. Transfer function W' 1 (z) exhibits some of the formant 

15 structure of the input speech signal. Thus, the masking property of the 
human ear is exploited by shaping the quantization error so that it has more 
energy in the formant regions where it will be masked by the strong signal 
energy present in these regions. The amount of weighting is controlled by 

the factors y } and y 2 - 

20 

The above traditional perceptual weighting filter 105 works well with 
telephone band signals. However, it was found that this traditional 
perceptual weighting filter 105 is not suitable for efficient perceptual 
weighting of wideband signals. It was also found that the traditional 
25 perceptual weighting filter 105 has inherent limitations in modelling the 
formant structure and the required spectral tilt concurrently. The spectral tilt 
is more pronounced in wideband signals due to the wide dynamic range 
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between low and high frequencies. The prior art has suggested to add a tilt 
filter into W{z) in order to control the tilt and formant weighting of the 
wideband input signal separately. 

A novel solution to this problem is, in accordance with the present 
5 invention, to introduce the preemphasis filter 1 03 at the input, compute the 
LP filter A(z) based on the preemphasized speech s{n), and use a modified 
filter W(z) by fixing its denominator. 

LP analysis is performed in module 104 on the preemphasized signal 
1 0 s(n) to obtain the LP filter A{z). Also, a new perceptual weighting filter 1 05 
with fixed denominator is used. An example of transfer function for the 
perceptual weighting filter 104 is given by the following relation: 

15 W{z) = a (z/ Y] ) / (i -y 2 z - 1 ) where 0<v <v <l 



A higher order can be used at the denominator. This structure substantially 
decouples the formant weighting from the tilt. 

20 

Note that because A(z) is computed based on the preemphasized 
speech signal s(n), the tilt of the filter 1/A(zS( } ) is less pronounced compared 
to the case when A(z) is computed based on the original speech. Since 
deemphasis is performed at the decoder end using a filter having the transfer 
25 function: 

/>-'(z)=i/(i-^- 1 ), 
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the quantization error spectrum is shaped by a filter having a transfer 
function W'\z)P~\z). When y 2 is set equai to /i, which is typically the case, 
the spectrum of the quantization error is shaped by a filter whose transfer 
5 function is ?/A(z/Y/), with A{z) computed based on the preemphasized 

speech signal Subjective listening showed that this structure for achieving 
the error shaping by a combination of preemphasis and modified weighting 
filtering is very efficient for encoding wideband signals, in addition to the 
advantages of ease of fixed-point algorithmic implementation. 

10 

Pitch Analysis: 

In order to simplify the pitch analysis, an open-loop pitch lag T OL is 
15 first estimated in the open-loop pitch search module 106 using the weighted 
speech signal sjn). Then the closed-loop pitch analysis, which is performed 
in closed-loop pitch search module 107 on a subframe basis, is restricted 
around the open-loop pitch lag T OL which significantly reduces the search 
complexity of the LTP parameters Tand b (pitch lag and pitch gain). Open- 
20 loop pitch analysis is usually performed in module 106 once every 10 ms 
(two subframes) using techniques well known to those of ordinary skill in the 
art. 

The target vector x for LTP (Long Term Prediction) analysis is first 
25 computed. This is usually done by subtracting the zero-input response s 0 of 
weighted synthesis filter W(z)/4(z) from the weighted speech signal s w (n). 
This zero-input response s 0 is calculated by a zero-input response calculator 
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108. More specifically, the target vector x is calculated using the following 
relation: 



5 

where x is the A/-dimensional target vector, s w is the weighted speech 
vector in the subframe, and s 0 is the zero-input response of filter W(z)/A(z) 
which is the output of the combined filter W(z)/A(z) due to its initial states. 

10 The zero-input response calculator 108 is responsive to the quantized 
interpolated LP filter A(z) from the LP analysis, quantization and interpolation 
calculator 104 and to the initial states of the weighted synthesis filter 
W(z)/A(z) stored in memory module 1 1 1 to calculate the zero-input response 
s 0 (that part of the response due to the initial states as determined by setting 

1 5 the inputs equal to zero) of filter W(z)/A(z). This operation is well known to 
those of ordinary skill in the art and, accordingly, will not be further described. 

Of course, alternative but mathematically equivalent approaches can 
be used to compute the target vector x. 

20 

A A/-dimensional impulse response vector h of the weighted 
synthesis filter W(z)/A(z) is computed in the impulse response generator 109 
using the LP filter coefficients A(z) and A(z) from module 104. Again, this 
operation is well known to those of ordinary skill in the art and, accordingly, 
25 will not be further described in the present specification. 

The closed-loop pitch (or pitch codebook) parameters b, Tand jare 
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computed in the closed-loop pitch search module 107, which uses the target 
vector x, the impulse response vector h and the open-loop pitch lag T OL as 
inputs. Traditionally, the pitch prediction has been represented by a pitch 
filter having the following transfer function: 

5 

1 / (1 -bz T ) 

where b is the pitch gain and T is the pitch delay or lag. In this case, the 
10 pitch contribution to the excitation signal u(n) is given by ba(n-T), where the 
total excitation is given by 

u{n) = bu(n-T)+gc k (n) 

15 

with g being the innovative codebook gain and c k (n) the innovative 
codevector at index fc 

20 This representation has limitations if the pitch lag Tis shorter than the 

subframe length N. In another representation, the pitch contribution can be 
seen as a pitch codebook containing the past excitation signal. Generally, 
each vector in the pitch codebook is a shift-by-one version of the previous 
vector (discarding one sample and adding a new sample). For pitch lags 

25 T>N, the pitch codebook is equivalent to the filter structure 1-bz r ) , and 
a pitch codebook vector v^ri) at pitch lag Tis given by 
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v T (n) = u (n-T) , n=0,...,AM. 

For pitch lags 7 shorter than N, a vector v^ri) is built by repeating the 
5 available samples from the past excitation until the vector is completed (this 
is not equivalent to the filter structure). 

In recent encoders, a higher pitch resolution is used which 
significantly improves the quality of voiced sound segments. This is 
10 achieved by oversampling the past excitation signal using polyphase 
interpolation filters. In this case, the vector v^ri) usually corresponds to an 
interpolated version of the past excitation, with pitch lag T being a non- 
integer delay (e.g. 50.25). 

15 The pitch search consists of finding the best pitch lag T and gain b 

that minimize the mean squared weighted error E between the target vector 
x and the scaled filtered past excitation. Error E being expressed as: 

20 E = \\x-by T \\ 2 

where y T is the filtered pitch codebook vector at pitch lag T: 

25 

n 

y T (") = v t (") * h(n) = Z v t (i)h(n-i) , n=0,...,/V-1. 
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It can be shown that the error E is minimized by maximizing the search 
criterion 



C=- 



* y 7 



1 0 where t denotes vector transpose. 

in the preferred embodiment of the present invention, a 1/3 
subsampie pitch resolution is used, and the pitch (pitch codebook) search is 
composed of three stages. 

15 

In the first stage, an open-loop pitch lag T OL is estimated in open-loop 
pitch search module 106 in response to the weighted speech signal sjn). 
As indicated in the foregoing description, this open-loop pitch analysis is 
usually performed once every 10 ms (two subframes) using techniques well 
20 known to those of ordinary skill in the art. 

In the second stage, the search criterion C is searched in the closed- 
loop pitch search module 107 for integer pitch lags around the estimated 
open-loop pitch lag T OL (usually ±5), which significantly simplifies the search 
25 procedure. A simple procedure is used for updating the filtered codevector 
y T without the need to compute the convolution for every pitch lag. 
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Once an optimum integer pitch fag is found in the second stage, a 
third stage of the search (module 107) tests the fractions around that 
optimum integer pitch teg. 

When the pitch predictor is represented by a filter of the form 
5 1/(1-bz T ), which is a valid assumption for pitch lags T>N, the spectrum of the 
pitch filter exhibits a harmonic structure over the entire frequency range, with 
a harmonic frequency related to 1/7. In case of wideband signals, this 
structure is not very efficient since the harmonic structure in wideband 
signals does not cover the entire extended spectrum. The harmonic 
10 structure exists only up to a certain frequency, depending on the speech 
segment. Thus, in order to achieve efficient representation of the pitch 
contribution in voiced segments of wideband speech, the pitch prediction 
filter needs to have the flexibility of varying the amount of periodicity over the 
wideband spectrum. 

15 

A new method which achieves efficient modeling of the harmonic 
structure of the speech spectrum of wideband signals is disclosed in the 
present specification, whereby several forms of low pass filters are applied 
to the past excitation and the low pass filter with higher prediction gain is 
20 selected. 

When subsample pitch resolution is used, the low pass filters can be 
incorporated into the interpolation filters used to obtain the higher pitch 
resolution. In this case, the third stage of the pitch search, in which the 
25 fractions around the chosen integer pitch lag are tested, is repeated for the 
several interpolation filters having different low-pass characteristics and the 
fraction and filter index which maximize the search criterion C are selected. 
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A simpler approach is to complete the search in the three stages 
described above to determine the optimum fractional pitch lag using only one 
interpolation filter with a certain frequency response, and select the optimum 
low-pass filter shape at the end by applying the different predetermined low- 
pass filters to the chosen pitch codebook vector v T and select the low-pass 
5 filter which minimizes the pitch prediction error. This approach is discussed 
in detail below. 

Figure 3 illustrates a schematic block diagram of a preferred 
embodiment of the proposed approach. 

10 

In memory module 303, the past excitation signal u(n) t n<0, is stored. 
The pitch codebook search module 301 is responsive to the target vector x, 
to the open-loop pitch lag T OL and to the past excitation signal u(n), n<0, from 
memory module 303 to conduct a pitch codebook (pitch codebook) search 

15 minimizing the above-defined search criterion C. From the result of the 
search conducted in module 301, module 302 generates the optimum pitch 
codebook vector v T . Note that since a sub-sample pitch resolution is used 
(fractional pitch), the past excitation signal u(n) t /?<0, is interpolated and the 
pitch codebook vector vy corresponds to the interpolated past excitation 

20 signal In this preferred embodiment, the interpolation filter (in module 301 , 
but not shown) has a low-pass filter characteristic removing the frequency 
contents above 7000 Hz. 

In a preferred embodiment, K filter characteristics are used; these 
25 filter characteristics could be low-pass or band-pass filter characteristics. 
Once the optimum codevector v r is determined and supplied by the pitch 
codevector generator 302, K filtered versions of v r are computed 
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respectively using K different frequency shaping filters such as 305®, where 
j=1, 2, ... , K. These filtered versions are denoted v f ® , where j=1, 2, ... , K. 
The different vectors v f & are convolved in respective modules 304®, where 
j=0, 1, 2, ... , K, with the impulse response h to obtain the vectors y®, where 
j=0, 1, 2, ... , K. To calculate the mean squared pitch prediction error for 
5 each vector y®, the value y® is multiplied by the gain b by means of a 
corresponding amplifier 307® and the value by® is subtracted from the target 
vector x by means of a corresponding subtractor 308®. Selector 309 selects 
the frequency shaping filter 305® which minimizes the mean squared pitch 
prediction error 

10 

e W=\\ x - b W y 0)§2 t j=1,2,...,K 

15 To calculate the mean squared pitch prediction error e® for each value of y®, 
the value y°> is multiplied by the gain b by means of a corresponding 
amplifier 307® and the value b®y«> is subtracted from the target vector x by 
means of subtractors 308®. Each gain b® is calculated in a corresponging 
gain calculator 306® in association with the frequency shaping filter at index 

20 j, using the following relationship: 

b U) =x tyU)/i y U)f 

25 

In selector 309, the parameters b, T, and /are chosen based on v T or 
vf> which minimizes the mean squared pitch prediction error e. 
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Referring back to Figure 1 , the pitch codebook index T is encoded 
and transmitted to multiplexer 112. The pitch gain b is quantized and 
transmitted to multiplexer 11Z With this new approach, extra information is 
needed to encode the index j of the selected frequency shaping filter in 
multiplexer 112. For example, if three filters are used (pO, 1, 2, 3), then two 
5 bits are needed to represent this information. The filter index information j 
can also be encoded jointly with the pitch gain b. 

Innovative codebook search: 

10 

Once the pitch, or LTP (Long Term Prediction) parameters b, T, and 
j are determined, the next step is to search for the optimum innovative 
excitation by means of search module 110 of Figure 1. First, the target 
vector x is updated by subtracting the LTP contribution: 

15 

x'=x-by T 



20 where b is the pitch gain and y T is the filtered pitch codebook vector (the 
past excitation at delay T filtered with the selected low pass filter and 
convolved with the inpulse response h as described with reference to Figure 

3). 

25 The search procedure in CELP is performed by finding the optimum 

excitation codevector and gain g which minimize the mean-squared error 
between the target vector and the scaled filtered codevector 
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£ = S x'- gHc k | 2 

where H is a lower triangular convolution matrix derived from the impulse 
5 response vector h. 

In the preferred embodiment of the present invention, the innovative 
codebook search is performed in module 110 by means of an algebraic 
codebook as described in US patents Nos: 5,444,816 (Adoul et al.) issued 
10 on August 22, 1995; 5,699,482 granted to Adoul et al., on December 17, 
1997; 5,754,976 granted to Adoul et al., on May 19, 1998; and 5,701,392 
(Adoul et al.) dated December 23, 1997. 

Once the optimum excitation codevector and its gain g are chosen 
15 by module 110, the codebook index k and gain g are encoded and 
transmitted to multiplexer 112. 

Referring to Figure 1, the parameters b, T t / A(z) t k and g are 
multiplexed through the multiplexer 112 before being transmitted through a 
20 communication channel. 

Memory update: 

25 In memory module 111 (Figure 1), the states of the weighted 

synthesis filter W(z)/A(z) are updated by filtering the excitation signal 
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u = gc k + bv T through the weighted synthesis filter. After this filtering, the 
states of the filter are memorized and used in the next subframe as initial 
states for computing the zero-input response in calculator module 108. 

As in the case of the target vector x, other alternative but 
5 mathematically equivalent approaches well known to those of ordinary skill 
in the art can be used to update the filter states. 

DECODER SIDE 

10 

The speech decoding device 200 of Figure 2 illustrates the various 
steps carried out between the digital input 222 (input stream to the 
demultiplexer 217) and the output sampled speech 223 (output of the adder 
221). 

15 

Demultiplexer 217 extracts the synthesis model parameters from the 
binary information received from a digital input channel. From each received 
binary frame, the extracted parameters are: 

20 - the short-term prediction parameters (STP) A(z) (once per frame); 

- the long-term prediction (LTP) parameters T, b f and ; (for each 
subframe); and 



25 



- the innovation codebook index k and gain g (for each subframe). 
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The current speech signal is synthesized based on these parameters as will 
be explained hereinbelow. 

The innovative codebook 218 is responsive to the index k to produce 
the innovation codevector which is scaled by the decoded gain factor g 
5 through an amplifier 224. In the preferred embodiment, an innovative 
codebook 218 as described in the above mentioned US patent numbers 
5,444,816; 5,699,482; 5,754,976; and 5,701,392 is used to represent the 
innovative codevector c k . 

10 The generated scaled codevector go k at the output of the amplifier 

224 is processed through a innovation filter 205. 

Periodicity enhancement: 

1 5 The generated scaled codevector at the output of the amplifier 224 

is processed through a frequency-dependent pitch enhancer 205. 

Enhancing the periodicity of the excitation signal u improves the 
quality in case of voiced segments. This was done in the past by filtering 

20 the innovation vector from the innovative codebook (fixed codebook) 218 
through a filter in the form 1/(1 -ebr 7 ) where £ is a factor below 0.5 which 
controls the amount of introduced periodicity. This approach is less 
efficient in case of wideband signals since it introduces periodicity over 
the entire spectrum. A new alternative approach, which is part of the 

25 present invention, is disclosed whereby periodicity enhancement is 
achieved by filtering the innovative codevector c k from the innovative 
(fixed) codebook through an innovation filter 205 (F(z)) whose frequency 
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response emphasizes the higher frequencies more than lower 
frequencies. The coefficients of F(z) are related to the amount of 
periodicity in the excitation signal u. 

Many methods known to those skilled in the art are available for 
5 obtaining valid periodicity coefficients. For example, the value of gain b 
provides an indication of periodicity. That is, if gain b is close to 1, the 
periodicity of the excitation signal u is high, and if gain b is less than 0.5, 
then periodicity is low. 

10 Another efficient way to derive the filter F(z) coefficients used in a 

preferred embodiment, is to relate them to the amount of pitch 
contribution in the total excitation signal u. This results in a frequency 
response depending on the subframe periodicity, where higher 
frequencies are more strongly emphasized (stronger overall slope) for 

15 higher pitch gains. Innovation filter 205 has the effect of lowering the 
energy of the innovative codevector c k at low frequencies when the 
excitation signal u is more periodic, which enhances the periodicity of the 
excitation signal u at lower frequencies more than higher frequencies. 
Suggested forms for innovation filter 205 are 

20 

(1) F(z)=1-oz~\ or (2) F(z) = -az+1-az _1 

25 where a or a are periodicity factors derived from the level of periodicity 
of the excitation signal u. 
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The second three-term form of F(z) is used in a preferred 
embodiment. The periodicity factor a is computed in the voicing factor 
generator 204. Several methods can be used to derive the periodicity 
factor a based on the periodicity of the excitation signal u. Two methods 
are presented below. 

5 

Method 1 : 

The ratio of pitch contribution to the total excitation signal u is first 
computed in voicing factor generator 204 by 

10 



15 



^ _ U V T V T _ n=0 



P 



U 'U 



lu 2 (n) 

n=0 



where v T is the pitch codebook vector, b is the pitch gain, and u is the 
excitation signal u given at the output of the adder 219 by 

20 u = gc k + bv T 

Note that the term bv T has its source in the pitch codebook (pitch 
codebook) 201 in response to the pitch lag T and the past value of u 
stored in memory 203. The pitch codevector v T from the pitch codebook 
25 201 is then processed through a low-pass filter 202 whose cut-off 
frequency is adjusted by means of the index /from the demultiplexer 217. 
The resulting codevector v T is then multiplied by the gain b from the 
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demultiplexer 217 through an amplifier 226 to obtain the signal bv T . 
The factor a is calculated in voicing factor generator 204 by 

<x ~ bounded by a < q 

5 

where q is a factor which controls the amount of enhancement (q is set 
to 0.25 in this preferred embodiment). 

Method 2: 

10 

Another method used in a preferred embodiment of the invention 
for calculating periodicity factor a is discussed below. 

First, a voicing factor r v is computed in voicing factor generator 204 

15 by 

20 where E v is the energy of the scaled pitch codevector bv T and E c is the 
energy of the scaled innovative codevector gc k . That is 

Ey = b 2 Vj t Vt = b 2 Yv* (n) 

n=0 

25 and 



W-1 

E r = g 2 cj c k = g 2 Tc' (n) 
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Note that the value of r v lies between -1 and 1 (1 corresponds to 
purely voiced signals and -1 corresponds to purely unvoiced signals). 

In this preferred embodiment, the factor a is then computed in 
voicing factor generator 204 by 

5 

a =0.125 (1 +r v ) 



10 which corresponds to a value of 0 for purely unvoiced signals and 0.25 
for purely voiced signals. 

In the first, two-term form of F(z), the periodicity factor o can be 
approximated by using o = 2a in methods 1 and 2 above, in such a 
15 case, the periodicity factor a is calculated as follows in method 1 above: 
o = 2qR p bounded by a < 2q. 

In method 2, the periodicity factor a is calculated as follows: 
o = 0.25 (1 + r v ). 

20 

The enhanced signal c f is therefore computed by filtering the 
scaled innovative codevector gc k through the innovation filter 205 (F(z)). 

The enhanced excitation signal u' is computed by the adder 220 

25 as: 
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u' - c f + fcvy 

Note that this process is not performed at the encoder 100. Thus, 
it is essential to update the content of the pitch codebook 201 using the 
5 excitation signal u without enhancement to keep synchronism between 
the encoder 100 and decoder 200. Therefore, the excitation signal u is 
used to update the memory 203 of the pitch codebook 201 and the 
enhanced excitation signal u' is used at the input of the LP synthesis filter 
206. 

10 

Synthesis and deemphasis 

The synthesized signal s' is computed by filtering the enhanced 
15 excitation signal u' through the LP synthesis filter 206 which has the form 
MA(z) 3 where A(z) is the interpolated LP filter in the current subfrarne. As 
can be seen in Figure 2, the quantized LP coefficients A(z) on fine 225 from 
demultiplexer 217 are supplied to the LP synthesis filter 206 to adjust the 
parameters of the LP synthesis filter 206 accordingly. The deemphasis filter 
20 207 is the inverse of the preemphasis filter 103 of Figure 1. The transfer 
function of the deemphasis filter 207 is given by 

D (z) = 1 / (1-mz~ 1 ) 

25 
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where fi is a preemphasis factor with a value located between 0 and 1 (a 
typical value is/z = 0.7). A higher-order filter could also be used. 

The vector s' is filtered through the deemphasis filter D(z) (module 
207) to obtain the vector s & which is passed through the high-pass filter 208 
5 to remove the unwanted frequencies below 50 Hz and further obtain s h . 

Oversampling and high-frequency regeneration 

0 The over-sampling module 209 conducts the inverse process of the 

down-sampling module 101 of Figure 1. In this preferred embodiment, 
oversampling converts from the 1 2.8 kHz sampling rate to the original 1 6 kHz 
sampling rate, using techniques well known to those of ordinary skill in the 
art. The oversampled synthesis signal is denoted s. Signal § is also 

5 referred to as the synthesized wideband intermediate signal. 

The oversampled synthesis § signal does not contain the higher 
frequency components which were lost by the downsampling process 
(module 101 of Figure 1) at the encoder 100. This gives a low-pass 
0 perception to the synthesized speech signal. To restore the full band of the 
original signal, a high frequency generation procedure is disclosed. This 
procedure is performed in modules 210 to 216, and adder 221 , and requires 
input from voicing factor generator 204 (Figure 2). 

5 In this new approach, the high frequency contents are generated by 

filling the upper part of the spectrum with a white noise properly scaled in the 
excitation domain, then converted to the speech domain, preferably by 
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shaping it with the same LP synthesis filter used for synthesizing the down- 
sampled signal § . 

The high frequency generation procedure in accordance with the 
present invention is described hereinbelow. 

5 

The random noise generator 213 generates a white noise sequence 
vf with a flat spectrum over the entire frequency bandwidth, using 
techniques well known to those of ordinary skill in the art. The generated 
sequence is of length AT which is the subframe length in the original domain. 
10 Note that N is the subframe iength in the down-sampled domain. In this 
preferred embodiment, A/=64 and A/-80 which correspond to 5 ms. 

The white noise sequence is properly scaled in the gain adjusting 
module 214. Gain adjustment comprises the following steps. First, the 
1 5 energy of the generated noise sequence vf is set equal to the energy of the 
enhanced excitation signal u* computed by an energy computing module 
210, and the resulting scaled noise sequence is given by 



A/-1 

E "' 2 ("> 

, /7=0,...,A/'-1. 

W'-1 

17=0 

25 The second step in the gain scaling is to take into account the high 

frequency contents of the synthesized signal at the output of the voicing 
factor generator 204 so as to reduce the energy of the generated noise in 



20 

w(n) = w'(n) 

\ 
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case of voiced segments (where less energy is present at high frequencies 
compared to unvoiced segments). In this preferred embodiment, measuring 
the high frequency contents is implemented by measuring the tilt of the 
synthesis signal through a spectral tilt calculator 212 and reducing the 
energy accordingly. Other measurements such as zero crossing 
5 measurements can equally be used. When the tilt is veiy strong, which 
corresponds to voiced segments, the noise energy is further reduced. The 
tilt factor is computed in module 212 as the first correlation coefficient of the 
synthesis signal s h and it is given by: 

10 

A/-1 

X) s h (n) s h (n -1) » conditioned by tilt ;> 0 and tilt > 

tm = JH2 

£ s„ 2 <»> 

15 

where voicing factor r v is given by 

20 

r v - (E„ - EJ I (E v + E c ) 

25 where E v is the energy of the scaled pitch codevector bv r and E c is the 
energy of the scaled innovative codevector gc h as described earlier. Voicing 
factor r v is most often less than tilt but this condition was introduced as a 
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precaution against high frequency tones where the tilt value is negative and 
the value of r v is high. Therefore, this condition reduces the noise energy for 
such tonal signals. 

The tilt value is 0 in case of flat spectrum and 1 in case of strongly 
5 voiced signals, and it is negative in case of unvoiced signals where more 
energy is present at high frequencies. 

Different methods can be used to derive the scaling factor g t from the 
amount of high frequency contents. In this invention, two methods are given 
1 0 based on the tilt of signal described above. 

Method 1 : 

The scaling factor g t is derived from the tilt by 

15 

g t = 1 - tilt bounded by 0.2 <; g t <;1.0 

For strongly voiced signal where the tilt approaches 1, g t is 0.2 and for 
20 strongly unvoiced signals g t becomes 1 .0. 

Method 2: 

The tilt factor g t is first restricted to be larger or equal to zero, then the 
25 scaling factor is derived from the tilt by 
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The scaled noise sequence w^produced in gain adjusting module 214 
is therefore given by: 

5 

w g = g t w. 

When the tilt is close to zero, the scaling factor g t is close to 1 , which 
does not result in energy reduction. When the tilt value is 1, the scaling 
1 0 factor g t results in a reduction of 1 2 dB in the energy of the generated noise. 

Once the noise is properly scaled (w g ), it is brought into the speech 
domain using the spectral shaper 215. In the preferred embodiment, this is 
achieved by filtering the noise w g through a bandwidth expanded version of 
1 5 the same LP synthesis filter used in the down-sampled domain (1/4(z/0.8)). 
The corresponding bandwidth expanded LP filter coefficients are calculated 
in spectral shaper 215. 

The filtered scaled noise sequence w f is then band-pass filtered to the 
20 required frequency range to be restored using the band-pass filter 216. In 
the preferred embodiment, the band-pass filter 216 restricts the noise 
sequence to the frequency range 5.6-7.2 kHz. The resulting band-pass 
filtered noise sequence z is added in adder 221 to the oversampled 
synthesized speech signal s to obtain the final reconstructed sound signal 
25 s out on the output 223. 
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Although the present invention has been described hereinabove by 
way of a preferred embodiment thereof, this embodiment can be modified at 
will, within the scope of the appended claims, without departing from the 
spirit and nature of the subject invention. Even though the preferred 
embodiment discusses the use of wideband speech signals, it will be 
obvious to those skilled in the art that the subject invention is also directed 
to other embodiments using wideband signals in general and that it is not 
necessarily limited to speech applications. 
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WHAT IS CLAIMED IS: 



1 . A pitch analysis device for producing an optimal set of pitch codebook 
parameters, comprising: 

a) at ieast two signal paths associated to respective sets of pitch 
5 codebook parameters, wherein: 

i) each signal path comprises a pitch prediction error 
calculating device for calculating a pitch prediction error of a pitch 
codevector from a pitch codebook search device; and 

ii) at least one of said two paths comprises a filter for 
10 filtering the pitch codevector before supplying said pitch 

codevector to the pitch prediction error calculating device of said 
one path; and 

b) a selector for comparing the pitch prediction errors calculated in said 
at least two signal paths, for choosing the signal path having the lowest 

15 calculated pitch prediction error, and for selecting the set of pitch codebook 
parameters associated to the choosen signal path. 

2. A pitch analysis device as defined in claim 1, wherein one of said at least 
two paths comprises no filter for filtering the pitch codevector before 

20 supplying said pitch codevector to the pitch prediction error calculating 
device. 

3. A pitch analysis device as defined in claim 1, wherein said signal paths 
comprises a plurality of signal paths each provided with a filter for filtering the 

25 pitch codevector before supplying said pitch codevector to the pitch 
prediction error calculating device of the same path. 
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4. A pitch analysis device as defined in claim 3, wherein the filters of said 
plurality of paths are selected from the group consisting of low-pass and 
band-pass filters, and wherein said filters have different frequency 
responses. 

5 5. A pitch analysis device as defined in claim 1, wherein each pitch 
prediction error calculating device comprises: 

a) a convolution unit for convolving the pitch codevector with a weighted 
synthesis filter impulse response signal and therefore calculating a convolved 
pitch codevector; 

10 b) a pitch gain calculator for calculating a pitch gain in response to the 

convolved pitch codevector and a pitch search target vector; 

c) an amplifier for multiplying the convolved pitch codevector by the pitch 
gain to thereby produce an amplified convolved pitch codevector; and 

d) a combiner circuit for combining the amplified convolved pitch 
1 5 codevector with the pitch search target vector to thereby produce the pitch 

prediction error. 

6. A pitch analysis device as defined in claim 5, wherein said pitch gain 
calculator comprises a means for calculating said pitch gain b U) using the 
20 relation: 

b Ul= x *y W/||y «J2 whefe j = a ^ ^ ^ K an(j K corresponds tQ a 

number of signal paths, 

25 and where x is said pitch search target vector, and y (J} is said convolved 
pitch codevector. 
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7. A pitch analysis device as defined in claim 1, wherein said pitch 
prediction error calculating device of each signal path comprises means for 
calculating an energy of the corresponding pitch prediction error, and 
wherein said selector comprises means for comparing the energies of said 
pitch prediction errors of the different signal paths and for choosing as the 

5 signal path having the lowest calculated pitch prediction error the signal path 
having the lowest calculated energy of the pitch prediction error. 

8. A pitch analysis device as defined in claim 5, wherein: 

a) each of said filters of the plurality of signal paths is identified by a filter 
10 index; 

b) said pitch codevector is identified by a pitch codebook index; and 

c) said pitch codebook parameters comprise the filter index, the pitch 
codebook index and the pitch gain. 



15 9. A pitch analysis device as defined in claim 1, wherein said filter is 
integrated in an interpolation filter of said pitch codebook search device, said 
interpolation filter being used to produce a sub-sample version of said pitch 
codevector. 



20 10. A pitch analysis method for producing an optimal set of pitch codebook 

parameters, comprising: 

a) in at least two signal paths associated to respective sets of pitch 

codebook parameters, calculating, for each signal path, a pitch prediction 

error of a pitch codevector from a pitch codebook search device; 
25 b) in at least one of said two signal paths, filtering the pitch codevector 

before supplying said pitch codevector for calculation of said pitch prediction 

error of said one path; and 
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c) comparing the pitch prediction errors calculated in said at least two 
signal paths, choosing the signal path having the lowest calculated pitch 
prediction error, and selecting the set of pitch codebook parameters 
associated to the choosen signal path. 

5 1 1 A pitch analysis method as defined in claim 10, wherein, in one of said 
at least two paths, no filtering of the pitch codevector is performed before 
supplying said pitch codevector to the pitch prediction error calculating 
device. 

10 12. A pitch analysis method as defined in claim 10, wherein said signal 
paths comprises a plurality of signal paths and wherein filtering the pitch 
codevector is performed in each of said plurality of signal paths before 
supplying said pitch codevector to the pitch prediction error calculating 
device of the same path. 

15 

13. A pitch analysis method as defined in claim 12, further comprising 
selecting the filters of said plurality of paths from the group consisting of low- 
pass and band-pass filters, and wherein said filters have different frequency 
responses. 

20 

14. A pitch analysis method as defined in claim 10, wherein calculating a 
pitch prediction error in each signal path comprises: 

a) convolving the pitch codevector with a weighted synthesis filter 
impulse response signal and therefore calculating a convolved pitch 

25 codevector; 

b) calculating a pitch gain in response to the convolved pitch codevector 
and a pitch search target vector; 
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c) multiplying the convolved pitch codevector by the pitch gain to thereby 
produce an amplified convolved pitch codevector; and 

d) combining the amplified convolved pitch codevector with the pitch 
search target vector to thereby produce the pitch prediction error. 

5 15. A pitch analysis method as defined in claim 14, wherein said pitch gain 
calculation comprises calculating said pitch gain b U) using the relation: 

b 01 -x V 01 /!!/ 01 !! 2 where j = 0, 1, 2, K, and K corresponds to a 

number of signal paths, 

10 

and where x is said pitch search target vector, and y* J) is said convolved 
pitch codevector. 

16. A pitch analysis method as defined in claim 10, wherein calculating said 
15 pitch prediction error, in each signal path, comprises calculating an energy 

of the corresponding pitch prediction error, and wherein comparing the pitch 
prediction error comprises comparing the energies of said pitch prediction 
errors of the different signal paths and choosing as the signal path having 
the lowest calculated pitch prediction error the signal path having the lowest 
20 calculated energy of the pitch prediction error. 

17. A pitch analysis method as defined in claim 14, wherein: 

a) identifying each of said filters of the plurality of signal paths by a filter 
index; 

25 b) identifying said pitch codevector by a pitch codebook index; and 

c) said pitch codebook parameters comprise the filter index, the pitch 
codebook index and the pitch gain. 
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18. A pitch analysis method as defined In claim 10, wherein said filtering the 
pitch codevector is integrated in an interpolation filter of said pitch codebook 
search device, said interpolation filter being used to produce a sub-sample 
version of said pitch codevector. 

5 1 9. An encoder having a pitch analysis device as in claim 1 for encoding a 
wideband input signal, said encoder comprising: 

a) a linear prediction synthesis filter calculator responsive to the 
wideband signal for producing linear prediction synthesis filter coefficients; 

b) a perceptual weighting filter, responsive to the wideband signal 
10 and the linear prediction synthesis filter coefficients, for producing a 

perceptually weighted signal; 

c) an impulse response generator responsive to said linear 
prediction synthesis filter coefficients for producing a weighted synthesis filter 
impulse response signal; 

15 d) a pitch search unit for producing pitch codebook parameters, 

said pitch search unit comprising: 

i) said pitch codebook search device responsive to the 
perceptually weighted signal and the linear prediction synthesis 
fitter coefficients for producing the pitch codevector and an 

20 innovative search target vector; and 

ii) said pitch analysis device responsive to the pitch 
codevector for selecting, from said sets of pitch codebook 
parameters, the set of pitch codebook parameters associated to 
the path having the lowest calculated pitch prediction error; 

25 d) an innovative codebook search device, responsive to the 

weighted synthesis fitter impulse response signal, and the innovative search 
target vector, for producing innovative codebook parameters; and 
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e) a signal forming device for producing an encoded wideband 
signal comprising the set of pitch codebook parameters associated to the 
path having the lowest pitch prediction error, said innovative codebook 
parameters, and said linear prediction synthesis filter coefficients. 

5 20. An encoder as defined in claim 19, wherein one of said at least two 
paths comprises no filter for filtering the pitch codevector before supplying 
said pitch codevector to the pitch prediction error calculating device. 

21. An encoder as defined in claim 19, wherein said signal paths comprises 
10 a plurality of signal paths each provided with a filter for filtering the pitch 

codevector before supplying said pitch codevector to the pitch prediction 
error calculating device of the same path. 

22. An encoder as defined in claim 21 , wherein the filters of said plurality of 
15 paths are selected from the group consisting of low-pass and band-pass 

filters, and wherein said filters have different frequency responses. 

23. An encoder as defined in claim 19, wherein each pitch prediction error 
calculating device comprises; 

20 a) a convolution unit for convolving the pitch codevector with the 

weighted synthesis filter impulse response signal and therefore calculating 

a convolved pitch codevector; 

b) a pitch gain calculator for calculating a pitch gain in response to the 

convolved pitch codevector and the pitch search target vector; 
25 c) an amplifier for multiplying the convolved pitch codevector by the pitch 

gain to thereby produce an amplified convolved pitch codevector; and 



WO 00/25298 



PCT/CA99/01008 



51 



d) a combiner circuit for combining the amplified convolved pitch 
codevector with the pitch search target vector to thereby produce the pitch 
prediction error. 

24. An encoder as defined in claim 23, wherein said pitch gain calculator 
5 comprises a means for calculating said pitch gain b (J} using the relation: 

b V^ily ^il 2 where j = 0, 1, 2, K, and K corresponds to a 

number of signal paths, 

10 and where x is said pitch search target vector, and y/ J) is said convolved 
pitch codevector. 

25. An encoder as defined in claim 19, wherein said pitch prediction error 
calculating device of each signal path comprises means for calculating an 

1 5 energy of the corresponding pitch prediction error, and wherein said selector 
comprises means for comparing the energies of said pitch prediction errors 
of the different signal paths and for choosing as the signal path having the 
lowest calculated pitch prediction error the signal path having the lowest 
calculated energy of the pitch prediction error. 

20 

26. An encoder as defined in claim 23, wherein: 

a) each of said filters of the plurality of signal paths is identified by a filter 
index; 

b) said pitch codevector is identified by a pitch codebook index; and 
25 c) said pitch codebook parameters comprise the filter index, the pitch 

codebook index and the pitch gain. 
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27. A encoder as defined in claim 19, wherein said filter is integrated in an 
interpolation filter of said pitch codebook search device, said interpolation 
filter being used to produce a sub-sample version of said pitch codevector. 

28. A cellular communication system for servicing a large geographical area 
5 divided into a plurality of cells, comprising: 

a) mobile transmitter/receiver units; 

b) cellular base stations respectively situated in said cells; 

c) a control terminal for controlling communication between the cellular 
base stations; 

10 d) a bidirectional wireless communication sub-system between each 

mobile unit situated in one cell and the cellular base station of said one cell, 
said bidirectional wireless communication sub-system comprising, in both the 
mobile unit and the cellular base station: 

i) a transmitter including an encoder for encoding a 
15 wideband signal as recited in claim 19 and a transmission circuit 

for transmitting the encoded wideband signal; and 

it) a receiver including a receiving circuit for receiving a 
transmitted encoded wideband signal and a decoder for decoding 
the received encoded wideband signal. 

20 

29. A cellular communication system as defined in claim 28, wherein one of 
said at least two paths comprises no filter for filtering the pitch codevector 
before supplying said pitch codevector to the pitch prediction error calculating 
device. 

25 

30. A cellular communication system as defined in claim 28, wherein said 
signal paths comprises a plurality of signal paths each provided with a filter 



WO 00/25298 



PO7CA99/01008 



53 

for filtering the pitch codevector before supplying said pitch codevector to the 
pitch prediction error calculating device of the same path. 

31. A cellular communication system as defined in claim 30, wherein the 
filters of said plurality of paths are selected from the group consisting of low- 

5 pass and band-pass filters, and wherein said filters have different frequency 
responses. 

32. A cellular communication system as defined in claim 28, wherein each 
pitch prediction error calculating device comprises: 

10 a) a convolution unit for convolving the pitch codevector with the 

weighted synthesis filter impulse response signal and therefore calculating 

a convolved pitch codevector; 

b) a pitch gain calculator for calculating a pitch gain in response to the 

convolved pitch codevector and the pitch search target vector; 
15 c) an amplifier for multiplying the convolved pitch codevector by the 

pitch gain to thereby produce an amplified convolved pitch codevector; and 
d) a combiner circuit for combining the amplified convolved pitch 

codevector with the pitch search target vector to thereby produce the pitch 

prediction error. 

20 

33. A cellular communication system as defined in claim 32, wherein said 
pitch gain calculator comprises a means for calculating said pitch gain b U) 
using the relation: 

25 h w =x f y ^/Hy^ll 2 where j = 0, 1, 2, K f and K corresponds to a 

number of signal paths, 
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and where x is said pitch search target vector, and y (J} is said convolved 
pitch codevector. 

34. A cellular communication system as defined in claim 28, wherein said 
pitch prediction error calculating device of each signal path comprises 

5 means for calculating an energy of the corresponding pitch prediction error, 
and wherein said selector comprises means for comparing the energies of 
said pitch prediction errors of the different signal paths and for choosing as 
the signal path having the lowest calculated pitch prediction error the signal 
path having the lowest calculated energy of the pitch prediction error. 

10 

35. A cellular communication system as defined in claim 32, wherein: 

a) each of said filters of the plurality of signal paths is identified by a 
filter index; 

b) said pitch codevector is identified by a pitch codebook index; and 
15 c) said pitch codebook parameters comprise the filter index, the pitch 

codebook index and the pitch gain. 

36. A cellular communication system as defined in claim 28, wherein said 
filter is integrated in an interpolation filter of said pitch codebook search 

20 device, said interpolation filter being used to produce a sub-sample version 
of said pitch codevector. 

37. A cellular mobile transmitter/receiver unit comprising: 

a) a transmitter including an encoder for encoding a wideband 
25 signal as recited in claim 19 and a transmission circuit for transmitting the 
encoded wideband signal; and 
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b) a receiver including a receiving circuit for receiving a transmitted 
encoded wideband signal and a decoder for decoding the received encoded 
wideband signal. 

38. A cellular mobile transmitter/receiver unit as defined in claim 37, wherein 
5 one of said at least two paths comprises no filter for filtering the pitch 

codevector before supplying said pitch codevector to the pitch prediction 
error calculating device. 

39. A cellular mobile transmitter/receiver unit as defined in claim 37, wherein 
10 said signal paths comprises a plurality of signal paths each provided with a 

filter for filtering the pitch codevector before supplying said pitch codevector 
to the pitch prediction error calculating device of the same path. 

40. A cellular mobile transmitter/receiver unit as defined in claim 39, wherein 
1 5 the filters of said plurality of paths are selected from the group consisting of 

low-pass and band-pass filters, and wherein said filters have different 
frequency responses. 

41 . A cellular mobile transmitter/receiver unit as defined in claim 37, wherein 
20 each pitch prediction error calculating device comprises: 

a) a convolution unit for convolving the pitch codevector with the 
weighted synthesis filter impulse response signal and therefore calculating 
a convolved pitch codevector; 

b) a pitch gain calculator for calculating a pitch gain in response 
25 to the convolved pitch codevector and the pitch search target vector; 
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c) an amplifier for multiplying the convolved pitch codevector by 
the pitch gain to thereby produce an amplified convolved pitch codevector; 

and 

d) a combiner circuit for combining the amplified convolved pitch 
codevector with the pitch search target vector to thereby produce the pitch 

5 prediction error. 

42. A cellular mobile transmitter/receiver unit as defined in claim 41 , wherein 
said pitch gain calculator comprises a means for calculating said pitch gain 
bf n using the relation: 

10 

b w =x f y w /ly w \\ 2 where j = 0, 7, 2, K, and K corresponds to a 

number of signal paths, 

and where x is said pitch search target vector, and y° ; is said convolved 
15 pitch codevector. 

43. A cellular mobile transmitter/receiver unit as defined in claim 37, wherein 
said pitch prediction error calculating device of each signal path comprises 
means for calculating an energy of the corresponding pitch prediction error, 

20 and wherein said selector comprises means for comparing the energies of 
said pitch prediction errors of the different signal paths and for choosing as 
the signal path having the lowest calculated pitch prediction error the signal 
path having the lowest calculated energy of the pitch prediction error. 

25 44. A cellular mobile transmitter/receiver unit as defined in claim 41, 
wherein: 
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a) each of said filters of the plurality of signal paths is identified by 
a filter index; 

b) said pitch codevector is identified by a pitch codebook index; 

and 

c) said pitch codebook parameters comprise the filter index, the 
5 pitch codebook index and the pitch gain. 

45. A cellular mobile transmitter/receiver unit as defined in claim 37, wherein 
said filter is integrated in an interpolation fitter of said pitch codebook search 
device, said interpolation filter being used to produce a sub-sample version 

10 of said pitch codevector. 

46. A cellular network element comprising: 

a) a transmitter including an encoder for encoding a wideband 
signal as recited in claim 19 and a transmission circuit for transmitting the 

1 5 encoded wideband signal; and 

b) a receiver including a receiving circuit for receiving a 
transmitted encoded wideband signal and a decoder for decoding the 
received encoded wideband signal. 

20 47. A cellular network element as defined in claim 46, wherein one of said 
at least two paths comprises no filter for filtering the pitch codevector before 
supplying said pitch codevector to the pitch prediction error calculating 
device. 

25 48. A cellular network element as defined in claim 46, wherein said signal 
paths comprises a plurality of signal paths each provided with a filter for 
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filtering the pitch codevector before supplying said pitch codevector to the 
pitch prediction error calculating device of the same path. 

49. A cellular network element as defined in claim 48, wherein the filters of 
said plurality of paths are selected from the group consisting of low-pass and 

5 band-pass filters, and wherein said filters have different frequency 
responses. 

50. A cellular network element as defined in claim 46, wherein each pitch 
prediction error calculating device comprises: 

10 a) a convolution unit for convolving the pitch codevector with the 

weighted synthesis filter impulse response signal and therefore calculating 

a convolved pitch codevector; 

b) a pitch gain calculator for calculating a pitch gain in response 

to the convolved pitch codevector and the pitch search target vector; 
15 c) an amplifier for multiplying the convolved pitch codevector by 

the pitch gain to thereby produce an amplified convolved pitch codevector; 

and 

d) a combiner circuit for combining the amplified convolved pitch 
codevector with the pitch search target vector to thereby produce the pitch 
20 prediction error. 

51 . A cellular network element as defined in claim 50 t wherein said pitch 
gain calculator comprises a means for calculating said pitch gain b a) using 
the relation: 

25 

h W= x V w 'Hy U1 \\ 2 where j = 0, 1, 2, K, and K corresponds to a 

number of signal paths, 
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and where x is said pitch search target vector, and y* J) is said convolved 
pitch codevector. 

52. A cellular network element as defined in claim 46, wherein said pitch 
5 prediction error calculating device of each signal path comprises means for 

calculating an energy of the corresponding pitch prediction error, and 
wherein said selector comprises means for comparing the energies of said 
pitch prediction errors of the different signal paths and for choosing as the 
signal path having the lowest calculated pitch prediction error the signal path 
10 having the lowest calculated energy of the pitch prediction error. 

53. A cellular network element as defined in claim 50, wherein: 

a) each of said filters of the plurality of signal paths is identified by 
a fitter index; 

15 b) said pitch codevector is identified by a pitch codebook index; 

and 

c) said pitch codebook parameters comprise the filter index, the 
pitch codebook index and the pitch gain. 

20 54. A cellular network element as defined in claim 46, wherein said filter is 
integrated in an interpolation filter of said pitch codebook search device, said 
interpolation filter being used to produce a sub-sample version of said pitch 
codevector. 



25 



55. In a cellular communication system for servicing a large geographical 
area divided into a plurality of cells, comprising: mobile transmitter/receiver 
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units; cellular base stations, respectively situated in said cells; and control 
terminal for controlling communication between the cellular base stations: 
a bidirectional wireless communication sub-system between each 
mobile unit situated in one cell and the cellular base station of said one cell, 
said bidirectional wireless communication sub-system comprising, in both the 
5 mobile unit and the cellular base station: 

a) a transmitter including an encoder for encoding a wideband 
signal as recited in claim 19 and a transmission circuit for transmitting the 
encoded wideband signal; and 

b) a receiver including a receiving circuit for receiving a transmitted 
1 0 encoded wideband signal and a decoder for decoding the received encoded 

wideband signal, 

56. A bidirectional wireless communication sub-system as defined in claim 
55, wherein one of said at least two paths comprises no filter for filtering the 

15 pitch codevector before supplying said pitch codevector to the pitch 
prediction error calculating device. 

57. A bidirectional wireless communication sub-system as defined in claim 
55, wherein said signal paths comprises a plurality of signal paths each 

20 provided with a filter for filtering the pitch codevector before supplying said 
pitch codevector to the pitch prediction error calculating device of the same 
path. 

58. A bidirectional wireless communication sub-system as defined in claim 
25 57, wherein the filters of said plurality of paths are selected from the group 

consisting of low-pass and band-pass filters, and wherein said filters have 
different frequency responses. 
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59, A bidirectional wireless communication sub-system as defined in claim 
55, wherein each pitch prediction error calculating device comprises: 

a) a convolution unit for convolving the pitch codevector with the 
weighted synthesis filter impulse response signal and therefore calculating 
a convolved pitch codevector; 
5 b) a pitch gain calculator for calculating a pitch gain in response 

to the convolved pitch codevector and the pitch search target vector; 

c) an amplifier for multiplying the convolved pitch codevector by 
the pitch gain to thereby produce an amplified convolved pitch codevector; 
and 

10 d) a combiner circuit for combining the amplified convolved pitch 

codevector with the pitch search target vector to thereby produce the pitch 
prediction error. 



60. A bidirectional wireless communication sub-system as defined in claim 
1 5 59, wherein said pitch gain calculator comprises a means for calculating said 
pitch gain b (J) using the relation: 

b «=x V w /Iy W B 2 where j = 0, 1, 2, K, and K corresponds to a 

number of signal paths, 



20 



and where x is said pitch search target vector, and y (J) is said convolved 
pitch codevector. 



61. A bidirectional wireless communication sub-system as defined in claim 
55, wherein said pitch prediction error calculating device of each signal path 
comprises means for calculating an energy of the corresponding pitch 
prediction error, and wherein said selector comprises means for comparing 
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the energies of said pitch prediction errors of the different signal paths and 
for choosing as the signal path having the lowest calculated pitch prediction 
error the signal path having the lowest calculated energy of the pitch 
prediction error. 

5 62. A bidirectional wireless communication sub-system as defined in claim 
59, wherein: 

a) each of said filters of the plurality of signal paths is identified by 
a filter index; 

b) said pitch codevector is identified by a pitch codebook index; 

10 and 

c) said pitch codebook parameters comprise the filter index, the 
pitch codebook index and the pitch gain. 

63. A bidirectional wireless communication sub-system as defined in claim 
15 55, wherein said filter is integrated in an interpolation filter of said pitch 
codebook search device, said interpolation filter being used to produce a 
sub-sample version of said pitch codevector. 
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